IEEE  JOURNAL  OF  OCEANIC  ENGINEERING,  VOL.  32,  NO.  3,  JULY  2007 


609 


Vertical  Directionality  of  Midfrequency  Surface 
Noise  in  Downward-Refracting  Environments 

Cathy  Ann  Clark,  Member,  IEEE 


Abstract — The  vertical  directionality  of  ambient  noise  due  to  sur¬ 
face  agitation  for  frequencies  between  2  and  5  kHz  propagated 
to  a  subsurface  receiver  has  a  characteristic  shape,  knowledge  of 
which  may  enhance  shallow- water  operations.  In  general,  the  noise 
level  is  highest  at  upward-looking  angles  and  attenuated  at  down¬ 
ward-looking  angles  depending  on  the  nature  of  the  bottom.  In  en¬ 
vironments  with  a  negative  profile  gradient,  the  noise  level  is  also 
greatly  reduced  in  a  low-angle  shadow  zone  or  “notch”  at  angles 
around  horizontal.  This  paper  reviews  the  character  of  vertical 
noise  directionality  by  examining  two  measured  data  sets  and  con¬ 
sidering  the  underlying  physical  mechanisms  that  drive  the  form  of 
the  distribution.  A  discussion  of  the  implications  of  vertical  noise 
directionality  for  design  and  operation  of  receiving  sonar  systems 
is  presented.  In  particular,  the  effect  of  mainlobe  beamwidth  and 
sidelobe  suppression  are  considered  along  with  the  directionality  of 
the  noise  field.  Finally,  an  overview  of  the  derivation  of  a  vertical 
noise  model  based  on  the  integrated  mode  method  of  propagation 
prediction  is  followed  by  model  reproduction  of  measurements. 

Index  Terms — Acoustics,  ambient  noise,  surface  noise,  under¬ 
water  sound  propagation,  vertical  noise  distribution. 


I.  Introduction 

IN  the  absence  of  near-shipping  interference,  the  undersea 
ambient  field  in  the  frequency  band  of  0.5-10  kHz  is 
dominated  by  noise  due  to  sea  surface  roughness  as  shown  in 
Fig.  1[1].  Although  near  and  distant  shipping  along  with  other 
undersea  acoustic  effects  such  as  internal  waves  are  also  factors 
in  the  ambient  noise  field,  this  paper  addresses  only  the  impact 
of  midfrequency  noise  resulting  from  sea  surface  agitation  near 
a  submerged  receiver  in  downward-refracting  environments. 
The  data  sets  studied  were  taken  in  scenarios  dominated  by  sea 
surface  noise,  and  the  shipping  components  of  the  noise  model 
described  in  Section  IV  are  not  discussed  herein. 

A  significant  amount  of  research  has  been  done  in  an  attempt 
to  quantify  and  explain  the  dependence  of  underwater  noise  on 
wind  speed  and  surface  agitation  [2]-[4],  yet  the  exact  nature  of 
the  coupling  of  energy  from  surface  agitation  due  to  wind  and 
sea  state  into  acoustic  pressure  waves  that  propagate  through  the 
undersea  environment  is  an  area  of  ongoing  research  [5]— [7] .  It 
is  believed  by  several  researchers  that  the  source  of  surface  noise 
is  breaking  waves  and  whitecaps  which  results  in  the  projection 
of  bubbles  into  the  water  column,  e.g.,  [4]— [7] . 

It  has  been  shown  that  dipole  and  monopole  models  of  the 
source  function  tend  to  apply  at  high  and  low  frequencies,  re- 


Manuscript  received  September  28,  2005;  revised  March  29,  2007;  accepted 
June  19,  2007. 

Associate  Editor:  D.  P.  Knobles. 

The  author  is  with  the  Naval  Undersea  Warfare  Center,  Newport,  Rl  02841 
USA  (e-mail:  clarkca@npt.nuwc.navy.mil). 

Digital  Object  Identifier  10. 11 09/ JOE. 2007. 903450 


140 


120 


100 


80 


40 


S. 


|  2° 


0 


Bubbles  and  Spray 


(Surface  Agitation) 


Agitation 


■aMM 

(Seismic  Background) 
—Turbulent-Pressure  Fluctuations 


(Surfnco  Wavos- 
— Second  Order  — 
Pressure  Effects) 


1  Precipitation 
hips,  Industrial,  Activity 


Biologies  — 

- ► 


Froquoncy  (Hz) 


10 


100  1,000  10,000  100,000 

.  Limb  of  Prevailing  Ncd«  - —  Earth  OuiVca  and  Erptosiw 

-  Wird-Depcndti*  Bubble  and  Spay  None  la/v-Fiequwcy  Very-StalSw-ttatefWhd 

—  —  Heavy  Preciptitco  Uiusl  Traffic  Nocs-Deep 

- Heavy  Traflc  Hod*  iiu»it«ii>-  iw^sk.iw. 


Fig.  1.  Acoustic  ambient  noise  in  the  ocean  (adapted  from  [1]  by  the  National 
Academy  of  Sciences). 


spectively  [8].  For  the  purposes  of  this  paper,  we  assume  that 
bubbles  from  the  surface  radiate  as  monopoles.  The  location  of 
a  monopole  source  below  a  pressure-release  surface  results  in  a 
dipole-like  propagation  pattern  from  the  perspective  of  a  subsur¬ 
face  receiver.  Propagation  of  each  source  through  the  environ¬ 
ment  and  superposition  of  pressure  (or  intensity)  at  the  position 
of  a  submerged  receiver  yield  the  noise  distribution  as  a  func¬ 
tion  of  vertical  angle  at  a  given  azimuth.  Further  discussion  of 
the  source  and  propagation  functions  is  given  in  Section  IV-A. 

The  ambient  noise  field  in  the  undersea  medium  as  perceived 
through  a  sonar  system  is  tactically  relevant  with  respect  to  a 
signal  arriving  simultaneously  at  the  same  sonar  receiver.  The 
3-D  convolution  of  the  noise  field  with  a  given  system  beam 
yields  the  beam  noise  level  associated  with  a  given  steer  angle. 
Repetition  over  a  set  of  vertically  directional  beams  yields 
the  distribution  of  beam  ambient  noise  as  a  function  of  beam 
depression/elevation  (D/E)  angle.  For  operational  use,  ambient 
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noise  from  the  surface  must  be  combined  with  other  effects, 
such  as  flow  noise,  which  depend  on  receiver  characteristics 
and  speed. 

This  paper  describes  two  sets  of  midfrequency  measurements 
of  vertical  noise  directionality  and  characterizes  the  nature  of 
the  observed  distributions.  The  basic  physical  acoustic  mech¬ 
anisms  determining  noise  directionality  are  reviewed  and  the 
implications  for  sonar  system  design  are  considered.  A  vertical 
noise  model  (V  Noise)  is  described  and  shown  to  reproduce  the 
shape  of  the  measured  data  sets  in  the  downward-refracting  en¬ 
vironments  studied. 

II.  Midfrequency  Data  Sets 

M  easurements  of  the  vertical  directionality  of  midfrequency 
surface  noise  have  been  made  in  a  handful  of  experiments.  The 
two  sample  distributions  included  herein  enable  the  general  na¬ 
ture  of  the  distribution  to  be  characterized  for  downward-re¬ 
fracting  environments,  i.e„  situations  in  which  sound  speed  de¬ 
creases  w  i  th  depth.  This  type  of  envi  ronment  was  chosen  i  n  each 
of  the  two  ex  peri  ments  to  obtai  n  measurements  of  the  I  ow-angl  e 
shadow  zone  or  “notch”  region  [9]. 

A.  Data  Set  1 — Gulf  of  Mexico,  June  1993 

The  first  of  the  two  data  sets  was  obtained  in  the  Gulf  of 
M  exico  (GOM )  in  J  une  1993  [10]  during  a  sea  test  sponsored  by 
the  Deployable  Acoustic  Sensor  System  (DASS)  project.  The 
site  of  the  test  was  about  100  nmi  west  of  Naples,  FL,  with 
the  receiver  moored  at  26°  19'  78"  N,  84°  23'  08"  W.  The 
measuring  array  was  part  of  the  U.S.  Naval  Research  Labora¬ 
tory  (NRL,  Washington,  DC)  digital  acquisition  buoy  system 
(NDABS)  and  was  provided  by  the  Stennis  Space  Center  Di¬ 
vision  of  the  Naval  Research  Laboratory  (NRL/SSC,  Stennis 
Space  Center,  MS)  [12],  It  was  a  31-element  vertical  line  array 
with  -28-dB  sidelobes  and  0.3-m  spacing,  and  thus,  a  design 
frequency  of  2500  Hz.  The  array  was  moored  with  its  center 
27  m  above  the  bottom  in  200  m  of  water  as  shown  in  Fig.  2. 


Individual  hydrophone  receptions  were  transmitted  to  an  instru¬ 
mented  pressure  vessel  where  they  were  digitized  at  a  rate  of 
6348  Hz  and  recorded  for  30  min  every  6  h  over  an  8-d  period. 

The  array  data  was  processed  by  performing  a  fast  Fourier 
transform  (F  FT )  on  data  from  each  hydrophone  with  a  transform 
sizeof  1024  and  a  Hanning  window,  resulting  in  a  3-dB  resolu¬ 
tion  of  approximately  9  Hz.  The  complex  data  were  then  beam- 
formed  in  each  FFT  bin  using  Taylor  shading  with  28-dB  side- 
lobes.  Finally,  the  beamformed  data  were  magnitude-squared 
and  block-averaged  over  128  transforms  (approximately  20  s). 

Sound-speed  profiles  (SSPs)  were  computed  from  sampled 
conductivity-temperature-depth  (CTD)  data  at  the  beginning 
and  end  of  the  test  and  a  fathometer  survey  indicated  a  slight 
bottom  slope  of  about  0.1°.  Information  provided  by  the 
Acoustic  Assessments  Section  of  NRL/SSC  indicated  that  the 
bottom  was  composed  of  muddy  sand  and  that  other  sources  of 
ambient  noise  such  as  shipping  activity  and  oil  drilling  were 
mi  ni  mal .  A  pi ot  of  the  S S P  at  the  beg i  nn i  ng  of  the  test  i s  show n 
in  Fig.  3  with  the  depth  of  the  receiving  array  center  indicated 
at  173  m. 

A  plot  of  beam  power  as  a  function  of  D/E  angle  (positive 
angles  are  up)  and  frequency,  measured  during  a  time  when 
there  was  no  ambient  contribution  from  biologies  or  shipping 
is  shown  in  Fig.  4.  The  low-angle  notch  is  evident  between  500 
and  2500  Hz  where  the  array  beamwidth  is  narrow  enough  to 
resolve  it.  Also  evident  in  Fig.  4  arestriation  patterns  which  are 
likely  due  to  a  resonance  effect  in  the  bottom  which  consisted  of 
a  layer  of  muddy  sand  over  hard  rock.  Frequencies  with  wave¬ 
lengths  which  would  be  trapped  in  the  sediment  layer  did  not  re¬ 
turn  to  the  water  column  for  multiples  of  some  base  frequency. 
A  slicethrough  Fig.  4  at  2473  Hz  isshown  in  Fig.  5.  This  single 
frequency  plot  of  beam  noise  versus  D/E  angle  will  be  used  in 
the  model  comparison  in  Section  V. 

B.  Measured  Data  Set  2  (DS2) — Tongue  of  the  Ocean,  Fall 
1988 

The  second  data  set  was  taken  from  an  experiment  in 
September/October  1988  in  the  Tongue  of  the  Ocean  (TOTO) 
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Fig.  4.  GOM  data  as  a  function  of  D/E  angle  and  frequency  (reproduced  from 
[10]). 


Fig.  6.  DAM  S  system  used  to  collect TOTO  data  (reproduced  from  [13]). 
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Fig.  5.  Slice  of  GOM  data  at  2473  Hz  (reproduced  from  [10]). 


in  the  Bahamas  at  a  depth  of  200  m  in  700  m  of  water  where 
the  bottom  depth  was  variable  in  the  near  proximity  of  the 
experiment  [8], 

The  measurement  system  was  a  deployable  acoustic-moni¬ 
toring  system  (DAM  S)  which  was  specifically  designed  to  mea¬ 
sure  the  vertical  distribution  of  sea  surface  sound.  The  wide¬ 
band  system  (40  to  4000  Hz)  was  composed  of  seven  nested 
arrays  of  seven  hydrophones  each,  with  -18-dB  sidelobes.  A 
sketch  indicating  the  hardware  components  of  the  system  is  re¬ 
produced  from  [13]  in  Fig.  6.  The  system  consists  of  “seven  oc- 
tavely  coalesced  four-wavelength  acoustic  line  antennas”  [13]. 
Each  antenna  section  transmitted  hydrophone  signals  to  a  sub¬ 
merged  data  logger  where  they  were  digitally  recorded.  A  full 
description  of  the  DA  M  S  system  and  subsequent  processing  of 
the  measured  noise  field  is  beyond  the  scope  of  this  paper.  A  n 


Fig.  7.  Illustration  of  the  location  of  the  hydrophone  array  relative  to  the  ocean 
surface  and  ocean  bottom  for  TOTO  data  (reproduced  from  [13]). 


Fig.  8.  TOTO  vertical  noise  as  a  dunction  of  D/E  angle  and  frequency  (repro¬ 
duced  from  [15]). 


overview  can  be  found  in  [13]  and  a  detailed  system  description 
in  [14].  The  vertical  noise  distribution  shown  in  Fig.  8  [15]  was 
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Fig.  10.  SSP  for  Norwegian  Sea  September  sample  environment. 


derived  from  the  ensemble  of  cross-spectral  density  functions 
for  all  possible  hydrophone  pairs  (see  [15]  and  [16]).  Nuttall 
[16]  has  assumed  that  the  inversion  method  used  captured  the 
complete  directional  spectrum.  Thus,  the  model  prediction  of 
the  vertical  directionality  of  the  4  kHz  noise  for  the  TOTO  data 
set  in  Section  V  assumes  omnidirectional  receive  sonar  beams. 

A  plot  of  the  receiver  orientation  with  respect  to  the  sloping 
bathymetry  in  the  TOTO  environment  is  reproduced  from  [13] 
in  Fig.  7,  and  the  SSP  is  shown  in  Fig.  9  with  the  depth  of 
the  receiving  array  indicated  at  200  m.  The  measured  vertical 
noise  distribution  as  a  function  of  frequency  and  vertical  angle 
is  shown  in  Fig.  8.  Due  to  isolation  of  the  TOTO  water  basin,  the 
ambient  noise  in  the  location  of  the  experiment  was  dominated 
by  local  sea  surface  conditions. 

III.  Sonar  Implications 

To  illustrate  the  propagation  effects  which  drive  the  subsur¬ 
face  vertical  directionality  of  surface  noise  and  to  investigate 
the  impact  of  beam  geometry  on  received  beam  noise,  a  loca¬ 
tion  in  the  Norwegian  Sea  with  a  receiver  at  a  depth  of  198  m 
in  a  water  column  of  1774  m  was  used  to  compute  the  vertical 
noise  distribution  and  beam  noise  for  a  sample  of  hypothetical 
beam  patterns.  The  downward- refracting  profile  for  the  sample 
environment  is  shown  in  Fig.  10(a)  with  an  enlarged  plot  of  the 
shallow  profile  in  Fig.  10(b).  The  bottom  was  modeled  using 
high-frequency  bottom-loss  (HFBL)  province  5  and  the  wind 
speed  was  13  kn,  the  historical  level  at  this  location  in  the  Nor¬ 
wegian  Sea  for  the  month  of  September. 

As  indicated  by  the  GOM  and  TOTO  data  sets  discussed 
in  Section  II,  surface  noise  at  a  submerged  receiver  in  down¬ 
ward- refracting  environments  is  vertically  nonisotropic,  with 
the  highest  level  at  upward-looking  beam  angles  and  a  lower 
level  at  downward-looking  beam  angles,  and  is  significantly  qui¬ 
eter  at  near-horizontal  angles.  These  three  regions  correspond, 
respectively,  to  noise  arriving  at  the  receiver  on  direct  paths 
from  the  surface,  noise  reflected  to  the  receiver  from  the  ocean 
bottom,  and  a  low-angle  notch  in  which  ray  paths  from  the  sur- 


Fig.  11.  Ray  traces  corresponding  to  angular  regions  of  a  sample  vertical  noise 
distribution. 


face  arriving  at  the  receiver  depth  will  originate  at  a  great  dis¬ 
tance  from  the  receiver,  and  hence,  be  significantly  attenuated. 
These  three  regions  and  their  effects  on  the  sample  vertical  noise 
curve  are  depicted  in  Fig.  11. 

The  sonar  array  and  beamformer  spatially  filter  the  incoming 
signal  and  noise  fields.  The  spatial  filters  or  beam  patterns  are 
formed  at  a  sufficient  number  of  vertical  and  azimuthal  steer  di¬ 
rections  to  sample  all  potential  signal  directions.  Both  thewidth 
of  thefiltermainlobe  and  the  degree  of  sidelobe  suppression  sig¬ 
nificantly  affect  the  received  beam  noise  level.  Fig.  12  depicts 
beam  noise  as  a  function  of  D/E  angle  for  a  set  of  beams  steered 
at  a  sample  of  vertical  D/E  angles  into  the  environmental  ver¬ 
tical  noise  distribution  as  shown.  The  calculation  of  total  beam 
noise  LN,  depicted  in  Fig.  12,  fora  single  steer  angle,  involves 
summing  the  dot  product  of  the  vertical  noise  level  as  a  func¬ 
tion  of  vertical  angle  V(6)  and  the  vertical  beam  response  as  a 
function  of  the  same  vertical  angles  B(8),  i.e., 

LN=  [  V(9)  ■  B(6)d0. 
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Fig.  12.  Beam  noise  level  versus  mainlobe  beamwidth  in  sample  environment. 
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Fig.  13.  Beam  noise  level  versus  sidelobe  level  in  sample  environment. 


The  vertical  axes  in  the  plots  in  Fig.  12(a)  and  (b)  are  vertical  di¬ 
rection  6,  while  the  vertical  axis  in  Fig.  12(c)  is  beam  steer  angle 
0B- The  beam  pattern  of  the  0°  D/E  angle  is  shown  with  1°,  2°, 
and  6°  beamwidths  plotted  in  red,  blue,  and  black,  respectively. 
The  received  beam  ambient  noise  levels  are  plotted  for  each  D/E 
angle  for  the  1°,  2°,  and  6°  mainlobe  beamwidths.  The  suppres¬ 
sion  of  sound  in  the  low-angle  notch  region  increases  with  the 
narrowing  of  the  main  beam.  The  comparative  noise  level  be¬ 
tween  upward-  and  downward-looking  angles  of  approximately 
15  dB,  however,  is  consistent  across  all  three  beamwidths.  These 
results  suggest  the  following:  1)  a  narrow  beam  relative  to  the 
width  of  the  notch  improves  notch  resolution  and  2)  the  dif¬ 
ference  in  gain  between  downward  and  upward  D/E  angles  is 
largely  independent  of  beamwidth. 

Since  the  beam  noise  level  at  a  given  D/E  angle  represents 
noise  arriving  at  all  vertical  angles,  sidelobe  suppression  also 
significantly  impacts  the  received  beam  noise  level.  Fig.  13 
shows  beam  noise  as  a  function  of  D/E  angle  for  a  set  of  beams 
against  the  same  hypothetical  noise  distribution  as  in  Fig.  12, 


but  with  sidelobe  suppressions  of  -13  and  -26  dB  shown 
in  blue  and  red,  respectively,  for  angles  between  ±40°.  As  in 
Fig.  12,  the  vertical  axes  in  Fig.  13(a)  and  (b)  are  vertical  direc¬ 
tion  0,  while  the  vertical  axis  in  Fig.  12(c)  is  beam  steer  angle 
6b ■  The  beam  noise  distributions  in  Fig.  13(c)  differ  in  that  the 
beam  with  greater  sidelobe  suppression  shows  a  lower  noise 
level  in  the  horizontal  notch  region  (over  10  dB  of  gain)  and  the 
gain  of  15  dB  in  downward-looking  over  upward-looking  D/E 
angles  is  consistent  across  beams  with  both  sidelobe  levels. 
These  results  suggest  that  increased  sidelobe  suppression  aids 
in  resolution  of  the  low-angle  notch. 

To  further  illustrate  the  effect  of  sidelobe  suppression  on 
beam  noise  directionality,  the  noise  arriving  on  the  mainlobe 
of  each  beam  versus  that  received  on  si  delobes  is  shown  in 
Fig.  14  for  a  beam  with  mainlobe  beamwidth  of  2°.  At  steeper 
angles  in  both  upward-  and  downward-looking  steer  angles, 
the  noise  curve  is  dominated  by  acoustic  energy  received  in  the 
mainlobe,  as  the  beams  in  these  regions  are  pointed  towards 
the  surface  or  bottom  where  the  noise  levels  are  highest.  For 


614 


IEEE  JOURNAL  OF  OCEANIC  ENGINEERING,  VOL.  32,  NO.  3,  JULY  2007 


D/E  Angle  (°) 


Fig.  14.  Beam  noise  received  on  mainlobe  and  sidelobes  in  sample 
environment. 


beams  steered  towards  the  horizontal  notch  region,  the  beam 
noise  is  dominated  by  energy  impinging  on  the  sidelobes  which 
are  directed  towards  angular  regions  of  higher  noise  level,  i.e., 
towards  the  surface  and  bottom.  The  change  in  sidelobe  level 
has  no  effect  on  the  upward-looking  D/E  angles  where  the  high 
level  of  noise  from  the  surface  is  received  on  the  mainlobe  of 
the  beam,  a  major  impact  in  the  notch  where  energy  is  highest 
in  the  sidelobes,  and  a  lesser  effect  on  downward-looking  D/E 
angles  where  the  highest  level  of  energy  is  again  received  on 
the  mainlobe,  but  has  been  attenuated  due  to  bottom  reflection. 

The  previous  discussion  deals  only  with  the  reception  of  am¬ 
bient  noise  by  the  sonar  system  beams.  In  addition,  there  are 
other  effects,  notably  flow  noise,  to  be  considered.  In  angular 
regions  where  the  ambient  noise  level  islow,  e.g.,  in  the  notch  re¬ 
gion  or  at  downward-looking  angles  over  an  absorbing  bottom, 
the  total  beam  noise  may  be  dominated  by  flow  noise.  Further¬ 
more,  the  tactical  relevanceof  vertical  directionality  of  ambient 
noise  depends  on  the  extent  to  which  it  masks  a  signal  of  in¬ 
terest.  Beam  noise  gain  associated  with  a  narrow  mainlobe  may 
be  accompanied  by  loss  of  signal  due  to  energy  splitting,  with  a 
resultant  decrease  in  signal-to-noise  ratio  (SNR).  The  best  SN  R 
may  be  in  a  direction  that  corresponds  to  neither  the  strongest 
signal  nor  the  quietest  noise. 

IV.  M  odel  Description 

The  noise  field  at  a  subsurface  receiver  due  to  wind  noise 
at  the  surface  is  modeled  by  integrating  the  effects  of  sound 
propagating  to  a  subsurface  receiver  from  multiple  monopole  or 
dipole  sources  distributed  throughout  an  area  of  the  sea  surface. 
This  approach  has  been  used  successfully  at  high  frequencies 
[17].  In  this  section,  we  discuss  the  theoretical  development  of 
a  vertical  noise  model  based  on  the  integrated  mode  method  of 
computing  propagation  documented  in  [18]. 

A.  Dipole  and  Monopole  Source  Factors 

In  ambient  noise  modeling,  the  source  of  surface  noise  is 
generally  treated  as  a  field  of  uniformly  distributed  subsurface 


Incremental  Surface  Area  dA-2  mdr 


Cross-sectional 
area  corresponding  to 
received  ray  bundle 


Ocean  bottom 

Fig.  15.  Overview  of  modeling  approach. 


Receiver 

Fig.  16.  Source  propagation  pattern. 


monopoles  and/or  dipoles  distributed  throughout  an  area  of  the 
sea  surface  [8],  [17],  as  depicted  in  Fig.  15.  The  model  com¬ 
putes  an  “intensity  density”  function  that  represents  power  per 
unit  of  cross-sectional  area  in  units  of  decibels  per  micropascal 
per  hertz  per  steradian.  Since  steradian  measure  corresponds  to 
the  surface  area  intersected  on  a  unit  sphere  with  a  total  surface 
area  of  An,  the  conversion  from  the  more  familiar  omnidirec¬ 
tional  decibel  level  to  decibels  per  steradian  is  -I01og10(47r) 
or  about  -11  dB. 

For  purposes  of  this  paper,  we  assume  a  monopole  source  at 
depth  z0  below  a  pressure- release  surface  and  a  reflected  point 
source  equidistant  above  the  surface  as  depicted  in  Fig.  16.The 
combination  of  direct  and  surface-  reflected  paths  results  in  a 
dipole-like  or  “doublet”  propagation  pattern  from  the  perspec¬ 
tive  of  a  subsurface  receiver. 

Assuming  plane-wave  propagation,  the  direct  path  contribu¬ 
tion  from  the  source  at  z0  is  given  by  a  plane  wave 

e~ikTr  e~ikz(zr-z0) 

where  kr  and  kz  are  the  horizontal  and  vertical  components  of 
the  wave  number  and  zr  is  the  depth  of  the  receiver.  The  re¬ 
flected  path  contribution  is  another  plane  wave 

■^e~ikrr  e-ik,(zr+zo) 

where  R  is  the  complex  surface  reflection  coefficient 
R  =  \R\e~iv. 

The  total  field  if)  of  the  two  waves  is 

if)  =  e~ikrr e-ikz(zr-z0)  _|_  e~ikz (z’'+z°) 

or 

if)  =  e-ikrre-ikzzr^eikzz0 •  +  Re-«’**»). 
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Fig.  18.  Cross-sectional  area. 


The  field  intensity  is  |^|2  or 

=  {e^z°  +Re~ik‘Zo)(erik^  +R*eik*m) 

where  the  *  denotes  complex  conjugation,  and  kr  and  kz  are 
assumed  to  be  real  quantities.  Simplifying,  the  field  intensity  is 

\ip\2  =  1  +  |i?.|2  -  2|i?|  cos(2kzz0)-  (1) 

If  \R\  =  1,  then  in  terms  of  angles  at  the  surface,  we  can  write 
kz  =  k0 sin 80,  and  the  intensity  reduces  to  the  dipole  pattern 

\ip\2  =  4sm2(kozosm0o).  (2) 

In  the  noise  model,  we  assume  azimuthal  invariance  and  use  a 
source  beam  B0  given  by 

Bo(60)  =1  +  |i?(#o)|2  -  2|i?,(#o)|  cos(2fc0^o  sin 8q).  (3) 


To  compute  the  cross-sectional  area  that  is  covered  by  the 
same  ray  bundle  at  range  r,  consider  the  geometry  of  Fig.  18. 
The  cross-sectional  area  is  Ar  =  2-Krsm8rdr.  The  amount  of 
ray  bundle  spreading  as  compared  to  that  on  the  unit  sphere  is 
given  by  the  ratio  Ar/A0.  The  sound  intensity,  which  expresses 
the  amount  of  spreading  of  the  ray  bundle  after  reaching  range 
r,  is  proportional  to  the  reciprocal  of  the  area  ratio.  Thus 

l  =  Ao 
Iq  Ar 

By  taking  the  reference  intensity  I0  to  be  unity,  then 

J.  COS  #0 

rsi"s''  (^) 

Since  we  express  angle  relationships  in  terms  of  6r,  we  use 
Snell’s  Law  to  replace  cos O0  by 


For  the  monopole  approximation,  the  two  paths  are  summed  c0  . 

randomly  in  phase  and  the  intensity  is  given  by  ^  cos  ' 


B0(80)  =  1  +  \R(80)\2.  (4) 

B.  Approximate  Ray  Theoretic  Analysis 

To  illustrate  the  method  in  a  simplified  context,  we  first  de¬ 
rive  an  expression  for  the  vertical  distribution  of  surface  noise 
at  a  submerged  receiver  based  on  ray  theoretic  considerations. 
U  nder  certain  conditions  (e.g.,  high-frequency  propagation  with 
8r  >  30°),  the  ray  theoretical  expression  is  valid,  but  for  small 
8r  (e.g.,  convergence  zone  propagation  with  8r  <  15°),  and  es¬ 
pecially  at  low  frequencies,  a  propagation  model  such  as  the 
integrated  mode  method  which  is  capable  of  computing  diffrac¬ 
tion  is  required. 

Consider  the  intensity  of  a  ray  bundle  of  width  dd0  emitted 
from  a  point  at  the  surface  at  an  angle  of  80  and  arriving  at  a 
subsurface  receiver  at  an  angle  of  8r  as  depicted  in  Figs.  17  and 
18. Theintensity  along  the  ray  bundleatagiven  horizontal  range 
is  related  to  the  number  of  steradians  per  unit  of  cross-sectional 
area  at  that  range,  where  the  steradi an  measure  of  the  ray  bundle 
i s  equal  to  the  surface  area  i ntersected  on  a  uni t  sphere.  The  area 
is  intersected  approximated  by  A0  =  2-k  cos  80d80  asillustrated 
in  Fig.  17. 


where  c0  and  cr  are  the  sound  speed  at  the  surface  and  receiver 
depth,  respectively.  Thus 

j Co  1  cos  8r 

cr  dr  sin  8r 

r - 

d80 

Now  consider  an  annular  region  on  the  surface  with  radius  r 
and  width  dr  as  shown  in  Fig.  15. 

Assuming,  for  the  moment,  a  density  of  noise  sources  with 
unit  strength  per  unit  of  area,  the  contribution  of  the  differential 
area  dA  at  the  receiver  is  IdA,  which  is  given  by  substituting 
dA  =  2%rdr  into  (5)  as 


IdA  = 


2ir 


Cq  COS  8 'f 

cr  sin  6r 


dOQ. 


Thus,  the  quantity  in  brackets  represents  an  “intensity  density” 
function  in  terms  of  6>0.  This  expression  can  be  extended  to  in¬ 
clude  azimuthal  variation  by  recognizing  that  the  2tt  factor  re¬ 
sults  from  integration  over  azimuthal  angle  cj>,  i.e., 

p2ir 

/  d(j)  =  27 r. 

Jo 
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TABLE  I 

Multipath/Angle  Correspondence 


Path 

Source  Angle 

Receive  Angle 

(1) 

down 

up 

(2) 

up 

up 

(3) 

down 

down 

(4) 

up 

down 

Thus,  the  “intensity  density”  V  in  both  horizontal  (0)  and  ver¬ 
tical  (</>)  directions  is  given  by 

V(8(h<j>)d60d<fr=  -^fdOod^. 

cr  sin  8r 

Since  we  wish  to  express  the  vertical  angle  density  in  terms  of 
the  receive  vertical  angle  8r,  we  perform  the  change  of  variable 

cos  8r  1 


to  obtain  density  in  8r  as 


This  result  corresponds  to  the  model  presented  in  [19]  and  [20], 


accounted  for,  the  coherent  combi  nation  of  these  two  paths  pro¬ 
duces  the  acoustic  field  of  a  dipole-like  or  doublet  source  at  the 
receive  point  as  discussed  in  Section  I V-A . 

In  the  V  Noise  model,  the  source  function  B0(80 )  is  included 
in  the  integrations  (7)  and  paths  /  =  2,4  are  omitted.  With  sur¬ 
face  and  bottom  reflection  coefficients  i?i  and  R2,  respectively, 
the  general  integral  form  (7)  with  kr) A(z0, kr)  re¬ 

placed  by  A  becomes 

hq  =  ~^=  f  2  B0(8)ARp1(R1R2y 
V2irr  Jki 

Xe-i(Kr+hTM(kr))dkr  (g) 

where  =  0 for 7  =  1  andp  =  -lforZ  =  3,  and  hr,iq  indicates 
vertical  phase  along  path  r  for  path  type  l  and  cycle  q. 

T  he  deri  vati  on  of  an  expressi  on  for  the  verti  cal  di  stri  buti  on  of 
surface  noise  at  the  receiver  proceeds  in  parallel  with  the  devel¬ 
opment  in  Section  IV-B,  i  .e.,  the  integrals  in  (8)  are  converted  to 
sound  intensity  and  the  contributions  of  individual  sources  are 
integrated  over  surface  area. 


D.  Calculation  of  Beam  Noise  and  Noise  Gain 

Noise  gain,  or  the  improvement  in  SN  R  due  to  reduction  in 
ambient  noise  due  to  the  directional  discrimination  of  a  given 
sonar  beam,  is  given  by  [21]  as  a  function  of  solid  angle  fl  as 


TVag 


101ogi0 


fN(Q)B(Q)dtt 


(9) 


C.  Integrated  M  ode  Formulation  of  Vertical  Noise  Distribution 

The  discussion  in  Section  IV-B  describes  the  sequence  of  cal¬ 
culations  used  to  construct  the  surface  noise  field  at  the  receiver 
usi  ng  si  mpl  e  ray  theory .  I  n  the  V  N  oi se  model  descri  bed  herei  n, 
the  same  development  applies  except  that  the  propagation  of  en¬ 
ergy  from  each  source  to  the  receiver  is  computed  by  the  inte¬ 
grated  mode  method  described  in  [18].  In  particular,  consider  an 
integral  corresponding  to  path  l  where  l  =  1,4  as  summarized 
in  Table  I  and  ray  path  cycle  q  where  q  =  0, 1, . . . 


y,:r  A  ( zr ,  )  ,4  ( ;_'o ,  kr ) 


xe-i(kTr+hr(kT))dkr 


where  B(Q)  is  the  receive  beam  pattern.  In  the  V  Noise  model, 
the  noise  distribution  is  normalized  so  that  integration  over  solid 
angle  S7  yields  the  midfrequency  omnidirectional  ambient  noise 
level  as  a  function  of  wind  speed  or  sea  state  and  frequency 
according  to  the  Wenz  curves  [1].  As  a  function  of  horizontal 
and  vertical  angles^  and  8,  the  total  omnidirectional  noise  level 
is  given  by 

L„  =  JJvm  cos  (8)d8d<t>  =  J  (10) 

Beam  ambient  noise  LEB  for  a  given  frequency  and  envi¬ 
ronment  is  computed  by  convolving  the  noise  V(8, </>)  with  the 
beam  pattern  B(8,  <j>)  at  each  steer  D/E  angle  0B  and  integrating 
the  result  over  azimuth  <f>,  i.e., 


where  r  is  a  horizontal  range,  kr  is  a  horizontal  wave  number, 
k\  and  k2  identify  upper  and  lower  turning  points  or  boundary 
reflections,  hr(kr )  is  an  elapsed  vertical  phase  along  path  r, 
A(zr, kr)  and  A(z0,kr)  are  wave  amplitudes  at  receiver  and 
source,  respectively,  and  a  ray  path  cycle  is  defined  as  propa¬ 
gation  from  one  upper  turning  point  or  reflection  to  the  next. 
For  convenience,  we  omit  boundary  losses  for  the  moment  so 
the  treatment  parallels  that  of  Section  IV-B. 

A  s  i  n  [18],  the  notation  for  path  type  is  such  that  l  =  1,2  rep¬ 
resent  propagation  paths  arriving  at  the  receiver  from  above  and 
l  =  3, 4  represent  propagation  paths  arriving  from  below.  Cycle 
q  =  O  represents  the  direct  path,  which  proceeds  from  a  given 
surface  source  at  depth  z0  to  the  receiver  at  depth  zr.  In  the  case 
depicted  in  Fig.  16,  the  combined  path  arrives  at  the  receiver 
from  the  upward  direction  and  thus  paths  l  =  1,2  correspond 
to  the  two  paths  in  Fig.  16.  With  a  suitable  boundary  condition 


LE b(8b)  =  I J  V(8,  <f>)  ■  B(0,  <f>)  cos(8)d8d<p.  (11) 

Finally,  total  beam  noise  LEBtot  is  the  power  sum  of  the 
ambient  noise  with  the  flow  noise  Le  for  a  given  receive  plat¬ 
form  and  speed.  The  SNR  for  a  given  target  is  computed  by 
comparing  the  noise  LE(8B)  to  the  received  signal  level  on  a 
beam-by-beam  basis. 

V.  Comparison  of  Model  Results 
TO  M  EASURED  DATA  SETS 

P I  ots  compari  ng  the  model  ed  and  measured  verti  cal  noisedis- 
tributions  for  the  GOM  and  TOTO  environments  of  Section  II 
are  shown  in  Figs.  19  and  20,  respectively.  In  both  plots,  the 
measured  values  are  shown  as  solid  dots  and  the  model  predic¬ 
tion  is  represented  by  a  solid  line. 
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Fig.  19.  M  odeled  versus  measured  noise  for  the  GOM  environment  (2473-Hz 
DASS  data). 


Fig.  20.  M  odeled  versus  measured  beam  noise  for  the  TOTO  environment. 


A  sea  state  level  of  2  was  assumed  for  the  GOM  data  set,  an 
estimate  obtained  from  [10],  and  the  modeled  beam  pattern  is 
shown  in  Fig.  21.  The  bottom  was  modeled  as  province  2  ac- 
cordi  ng  to  the  H  F  B  L  data  base  from  the  0  ceanographi c  and  A  t- 
mospheric  M  aster  Library  [22],  which  is  valid  for  frequencies 
between  1.5  and  4.0  kHz,  because  it  most  closely  matches  the 
bottom  loss  predicted  by  modified  acoustic  bottom  loss  evalua¬ 
tion  (ABLE)  [11]  province  1  used  in  [10],  An  unidentified  noise 
source  resulted  in  a  notch  floor  in  the  measurements  of  about 
25  dB  which  is  not  predicted  by  the  model. 

For  the  TOTO  data  set,  a  wind  speed  of  25  kn  and  HFBL 
province  2  were  used.  Since  the  model  normalizes  the  vertical 
noise  distribution  to  match  the  overall  omnilevel  predicted  for  a 
given  sea  state  and  frequency  according  to  the  Wenz  curves  [1], 
a  wind  speed  of  25  kn  and  a  bottom  type  of  HFBL  province  2 
resulted  in  a  match  to  the  level  of  the  measured  data.  Although 
there  is  no  basis  for  these  assumptions  other  than  matching  to 
the  data,  the  shape  of  the  noise  distribution  has  been  success¬ 
fully  reproduced  by  the  model.  The  model  eddistributi  on  isF((9) 


Fig.  21.  N  DABS  beam  pattern  at  steer  angle  0°  (horizontal). 


rather  than  LE (0B)  because  no  beam  was  appl  i  ed.  T  he  i  nversi  on 
method  used  in  processing  the  data  is  believed  to  have  recovered 
thefull  directional  spectrum  of  the  measured  noise  field  [16],  as 
discussed  in  Section  II. 

The  source  model  B0(9)  is  the  effective  dipole  (3)  of 
Section  IV-A,  which  represents  the  interference  pattern  be¬ 
tween  a  monopole  located  below  the  pressure- release  surface 
and  its  image  above  the  surface.  To  control  phase  interference 
patterns  as  frequency  is  increased,  the  product  of  horizontal 
wave  number  k  times  source  depth  zs,  which  appears  in  (3) 
for  the  dipole,  has  been  set  equal  to  a  constant  value  of  1/3, 
a  constant  previously  chosen  by  comparison  with  measured 
data.  To  test  the  effectiveness  of  this  assumption  with  respect 
to  the  current  data  sets,  the  predicti  ons  for  the  G  O  M  and  TOTO 
environments  were  repeated  using  a  number  of  possible  source 
functions.  The  six  functions  modeled  are  as  follows: 

1)  function  proposed  by  Becken  [23]  given  by 


where 


Bo(0o)=smx9o 


2,  90  <  20° 

fS,  0o>  20° 


(12) 


2)  dipole  pattern  suggested  by  Carey  [24]  given  by 

Bo(60)  =  1  +  |i?((9o)|2  -  2|i?,(#0)|  cos(2kzs sin #0) 


where  A;  is  a  horizontal  wave  number  and  the  source  depth 
zs  is  set  equal  to  A/4  where  A  is  a  wavelength,  i.e.,  the 
product  of  kzs  is  always  equal  to  7r/2; 

3)  sin(#0); 

4)  sin2(#0); 

5)  dipole  used  herein,  i.e., 

Bo(60)  =  1  +  |-R($o)|2  -  2|i?(#0)|  cos(2kzs sin 9q) 


where  zs  =  1  yard  and  kzs  =  1/3; 

6)  monopole  source  function  1  +  \R(60)\2- 
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Fig.  22.  Comparison  of  GOM  measurements  with  modeled  vertical  noise  for 
various  source  functions. 


Fig.  23.  Comparison  of  TOTO  measurements  with  modeled  vertical  noise  for 
various  source  functions. 


Although  the  source  function  used  currently  has  no  theoret¬ 
ical  basis  to  defend  it,  the  results  of  the  comparison  shown  in 
Figs.  22  and  23  for  GOM  and  TOTO,  respectively,  indicate  that 
the  best  agreement  with  the  two  data  sets  studied  is  achieved 
by  its  use.  Further  research  is  needed  to  derive  a  more  rigorous 
source  function. 

VI.  Summary 

In  summary,  it  is  noted  that  the  undersea  environment  can  pro¬ 
vide  a  strong  vertical  noise  directionality,  especially  when  the 
SSP  is  downward  refracting.  The  highest  noise  arrives  at  a  sub¬ 
surface  receiver  from  upward  directions.  A  significant  quiet  re¬ 
gion  exists  at  near-horizontal  directions  in  downward-refracting 
environments  and  the  noise  from  downward-looking  directions 
is  attenuated  according  to  the  nature  of  the  bottom.  The  major 
features  of  the  vertical  noise  distribution  have  been  discussed 
in  terms  of  the  physical  mechanisms  which  cause  the  direc¬ 
tionality.  The  impact  on  beam  geometry  for  a  subsurface  re¬ 
ceiver  with  respect  to  noise  suppression  and  signal  detection 


have  been  investigated.  The  VNoise  model  successfully  repro¬ 
duces  the  shape  of  the  measured  midfrequency  vertical  noise 
distributions  studied. 
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